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This  paper  considers  a  cantilever  beam  which  can  move  in  and  out 
of  a  rigid  frictionless  hole.  Hamilton's  principle  is  used  for 
the  development  of  two  governing  equations,  and  boundary  condi¬ 
tions  defined  over  a  moving  boundary.  A  transformation  of 
coordinates  is  used  to  transform  the  problem  into  a  system  of 
two  nonlinear  p.d.e.  defined  over  a  fixed  region.  A  discussion 
of  some  mathematical  and  physical  considerations  is  presented. 
Finally,  a  finite  element  formulat ion  of  the  problem  is  presente 


UNCLASSIFIED 


CONTENTS 


Page 


1.  Introduction  ------------------  1 

2.  Analysis  of  Equations  -------------  6 

3.  Method  of  Solution  ---------------12 

4.  Figures  ------  17 

REFERENCES - - 19 

INITIAL  DISTRIBUTION  LIST  -  20 


I  !  I 


VIBRATION  OF  A  CANTILEVER  BEAM  THAT 


SHOES  AXIALLY  IN  A  RIGID  FRICTIONLESS  HOLE 


INTRODUCTION 

Consider  a  cantilever  beam  that  fits  snugly  into  a  frictionless  hole, 
Fig.  I.  Axes  x,y  are  fixed.  A  force  F(t)  moves  the  beam  axially.  Hence, 
the  length  L  of  the  beam  outside  of  the  hole  is  a  function  of  time  t.  The 
lateral  deflection  of  the  beam  is  y(x,t)  ,  where  for  x  <  0,  y  •  0.  The 
axial  velocity  of  the  beam  is  ll  •  dL/dt. 

By  Hamilton's  principle, 

(5T  ♦  5W)  dt  •  0  (1) 

o 


where  T  denotes  the  kinetic  energy  of  the  system  and  W  denotes  the  work 
of  external  and  internal  forces.  The  kinetic  energy  is 


(2) 


where  the  subscript  t  denotes  partial  differentiation,  p  •  the  mass  of  the 

beam  per  unit  length  ■  constant,  and  L  -  the  length  of  the  beam  when  t  ■  0. 

o 

The  strain  energy  of  the  beam  (the  negative  of  the  work  of  internal  forces 
for  a  conservative  system)  is 


1 


(3) 


t 


o 


where  E  -  the  modulus  of  elasticity  of  the  beam  and  I  -  the  second  moment 
of  the  beam  cross-sectional  area.  Hence,  the  virtual  work  6W  is 


5U  -  F6L  -  5U 


(*> 


where  SU  ■  the  first  variation  of  the  strain  energy  (Eq.  3).  Because  of 
the  inertia,  the  axial  movement  of  the  beam  tends  to  bend  the  beam  by  beam- 
column  action.  This  effect  is  neglected  here,  but  it  could  be  included. 
Equations  (I),  (2),  (3),  and  (I*)  yield 


/ 


[pLoL3L  +  i  pyt2(l,t)  5L  ♦  FfiL  -  1  Ely^  (L,t)  Si 


I  ..  2 


L 

*  °  f  yt5ytdx  "  El  j 


dx  -  El  /  y  5y  dx]  dt  »  0 
'  xx  'xx 


(5) 


If  L(t)  is  considered  as  given,  the  terms  with  the  factors  6l  and  <5L  merely 
determine  F(t)  to  give  the  prescribed  motion  L(t).  The  Euler  equation  for 
the  lateral  motion  of  the  beam  is 


Ely  +  py  "0 

' xxxx  7tt 


(6) 


The  natural  boundary  conditions  are 


yxx(L*t}  "  yxxxa-c)  -  0 


(7) 
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T 


The  forced  boundary  conditions  are 


y (o,t)  -  yx(o,t)  •  0  (8) 

Adopting  Hamilton's  viewpoint,  we  may  suppose  that  y(x,o)  and  y(x,tj) 
are  prescribed  (Fig.  2).  The  problem  is  to  determine  y(x,t)  in  region  R 
(Fig.  2)  subject  to  the  differential  equation  (Eq.  6)  and  the  indicated 
boundary  conditions  (Fig.  2). 

The  curved  boundary  at  the  right  side  of  region  R  is  troublesome.  It 
is  possible  to  transform  coordinates  so  that  this  boundary  becomes  straight. 

Let 


^  ,  t  •  tU,q) 


(9) 


Regard  y  as  a  function  of  (5,n).  Then, 


y  ■  y(5,n),  5  -  S(x, t) ,  n  -  n(x,t) 


do) 


and 


■  i 

(~J'x  ■ 


on 


Here,  we  regard  y  as  being  a  function  of  U.n)  . 
Now  let 


s  L 


r  .  t  *  n 


0  2) 


3 


we  write  ,  etc.,  we  imply  y  -  y(g,n).  Hence,  |j  -  (||)  (I)  . 


end  s imi lar ly , 


and  since 


we  f i nd 


and  hence, 


,4  ,  .4 

3  y  _  I  9  y 

4  T  *  — IT 

3 xH  LH  3g* 


(45.)  -  -  £k  -  -  5k  (Hi) 

«  *  l2  L  *  '3t  ; 


;  .  -(ii\  sL  +  ( h 

v  L  '3n  g 


(ii)  «  -(iiil)  ik  .  rll\  4  32y 

Hn  ^-2}  L  y-iL+^ 

(li.)  „  .  32y  gL  ,9y.  gi.  ,9vv  gi' 
3n  g  3g3n  L  3g  n  T”  3g  n  *l" 


Thus,  by  the  third  of  Eqs.  (II),  Eqs.  (12),  Eqs.  (14).  and  Eqs.  (15),  we 


t*>* '  £  -  <$£■ *  £&  i-fr 


-;2  .2 


L  3n 


4 


and  since  n  ■  t,  we  obtain 


,2  ..2  .2:2  .  2  .2 

(M),  V2 

3tZ  3C  LZ  ^  L 


-  2lSk  ♦  i 


3"y  £L_ 

3£3t  L 


K  L 


3t 


(16) 


in  which  y  ■  y(£,t)  on  the  right  hand  side  of  Eq.  (16). 
Consequently,  Eqs.  (6),  (13)  and  (16)  yield 


..  J  ,2  .2-2 

El  3  v  *  -,[.Lx  i_L 

L'  ZC  3C2  L2 


*2fe  *£ 


2  ••  2 
,  3  y  £1  ik  * 

•  2  307  It  *  st  It  * 


(17) 


where  y  is  now  regarded  as  a  function  of  Z  and  t.  Equations  (7)  and  (8) 
become 


y..(l ,t)  -  y---(l , t)  -  0 

(18) 

y  (0,t)  -  y.  (0,t)  -  0 


Equation  (17)  is  much  more  complicated  than  Eq.  (6).  However,  the 
region  R  is  simpler  (Fig.  3)-  The  problem  remains  a  linear  boundary-value 
problem  if  in  Eq.  (17)  we  regard  l(t)  as  a  given  function  of  time  t.  If 
we  consider  L(t)  as  unknown,  we  must  obtain  another  equation  which  defines 
l(t).  This  equation  is  obtained  from  Eq.  (5)  from  the  terms  with  the 
factors  5l  and  5L.  Hence,  we  find 
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(19) 


Equation  (20)  is  subject  to  the  initial  conditions 
L(0)  -  Lo 

L  (0)  -  0 


(21) 


where  l  is  the  initial  length  of  the  beam  at  time  t  *  0,  (see  Fig.  1). 
o 

Thus,  if  y(x,t)  and  L(t)  are  considered  as  the  unknowns,  Eqs.  (17) 
and  (20)  subject  to  the  boundary  conditions  of  Eqs.  (18)  and  (21),  respec¬ 
tively,  represent  a  complicated  nonlinear  boundary  value  problem. 

ANALYSIS  OF  EQUATIONS 

(a)  L ( t )  Prescribed.  If  L(t)  is  given,  then  the  problem  reduces  to 
seeking  a  solution  of  Eq.  (17)  subject  to  the  boundary  conditions,  Eqs.  vl8). 
The  problem  is  still  of  considerable  difficulty.  However,  the  problem  is 
a  linear  boundary-val ue  problem.  The  nature  of  the  solutions  of  Eq .  (17) 
is  unknown.  Indeed,  one  must  consider  whether  or  not  a  unique  solution  of 
the  boundary-va I ue  problem  exists.  Also,  is  there  a  unique  solution  of 
the  corresponding  initial  value  problem  in  which  the  condition  y  ■  Fj(£) 
on  t  *  t|  is  discarded  and  the  functions  y(C.O)  and  y^U.O)  are  prescribed? 
One  way  to  approach  such  questions  is  through  the  calculus  of  finite 
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differences,  in  which  the  differencial  equation  is  replaced  by  a  linear 
partial  difference  equation.  The  existence  of  a  solution  of  the  difference 
equation  is  a  question  of  linear  algebra.  The  most  practical  way  of 
developing  numerical  solutions  of  Eqs..  (17)  and  (18)  is  an  open  question. 

One  possible  approach  is  through  the  calculus  of  finite  differences,  or  a 
combination  of  the  calculus  of  finite  differences  with  the  t<unge-K.utta 
method.  Another  approach  is  through  finite-element  methods. 

The  problem  is  simplified  greatly  if  the  axial  motion  is  so  slow 
that  time  derivatives  of  L(t)  are  negligible.  For  example,  if  the  beam  is 
vibrating  in  a  natural  mode  and  it  is  then  drawn  slowly  into  the  hole, 
does  it  continue  to  vibrate  in  a  natural  mode?  How  do  the  frequency  and 
amplitude  vary  with  l  under  these  conditions?  For  a  rigid,  frictionless 
hole,  the  total  mechanical  energy  is  constant,  since  no  energy  is  imparted 
at  the  root.  This  criterion  serves  to  determine  how  the  amplitude  varies. 
However,  we  should  not  have  to  apply  this  criterion,  since  constancy  of 
mechanical  energy  should  follow  automatically  from  the  solution  of  Eq.  (17). 

If  L  and  L  are  negligible,  Eq.  (17)  reduces  to 


34V  + 

3S4  £l  It1 


0 


(22) 


Equation  (22)  is  much  simpler  than  Eq.  (17).  However,  the  fact  that  L  is 
a  prescribed  function  of  t  still  complicates  it.  There  is  a  question  as 
to  whether  the  character  of  the  solution  is  altered  by  the  reduction  of 
Eq.  (17)  to  (22).  Again,  does  the  initial  value  problem  posed  by 
Eqs.  (18)  anc  (.22 )  have  a  unique  solution?  If  so,  the  constancy  of  mechanical 
energy,  as  a  separate  principle,  is  not  needed  to  determine  the  dependency 
of  amplitude  upon  L.  In  fact,  constancy  of  mechanical  energy  should  be  a 
deducible  consequence  of  Eqs.  (13)  and  (22). 
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Unlike  Eq.  (17),  Eq.  (22)  allows  a  theory  of  natural  modes.  For 
example,  consider  solutions  of  Eq.  (22)  of  the  form  y  =  f(£)g(t).  Then, 
Eq.  (22)  yields 


Hence , 


f""  OL^.  £1 

f  El  g 


0 


f""  ,4  Dlk  g"  04 

—  "  6  •  IT  S 


(23) 


in  which  p  is  a  constant.  Then, 


34EI. 


IT  9 


OL 


(24) 


Unfortunately,  Eq.  (24)  does  not  yield  simple  harmonic  motion,  because  L 
is  a  function  of  t.  However,  the  equation 


f  I  II  I 


-  S4f 


(25) 


can  be  integrated.  With  the  root  conditions  (Eqs.  18)  f(o)  =  f'(o)  =  0, 
Eq.  (25)  yields 

f(c|)  *  A(sinhbC  -  sinp£)  +  B(coshpi;  -  cospE)  (26) 


Now,  the  free-end  conditions  (Eqs.  18)  f"(l)  ■  f"*(l)  -  0  yield 


A(sinhS  +  sin6)  +  8 (coshS  +  cosp)  ■  0 
A(coshj3  +  cosSl  +  S(sinh8  -  sing)  ■  0 


and  the  vanishing  of  the  determinant  of  Eqs.  (27)  yields 


(27) 


3 


coshB  cosB  ■  -1 


(28) 


These  results  are  exactly  like  that  for  natural  modes  of  a  cantilever  with 
fixed  length.  The  roots  of  Eq.  (28)  are  constants  8, ,  S2»  Sj  ,  •  •  •  •  They 
have  been  tabulated  by  Young  and  Felgar  (1). 

As  is  customary,  let  B  *  1  and  A  »  _an’  Then,  by  Eq.  (27), 

coshB  +  cosB_  ,  % 

_ _ -2 - r-Jl  (29) 

n  sinhB  +  si nB_ 
n  n 

Consequently,  a] ,  a2>  •••  are  constants,  independent  of  l.  The  values  of 

a  have  been  tabulated  by  Young  and  Felgar.  Equation  (26)  gives  the  natural 
n 

modes , 

f  (£)  =  coshB  5  -  cosB  £  -  a  (sinhB  l  ~  sinB  E)  (30) 

n  n  n  n  n  n 

The  functions  f  (£)  are  orthogonal  in  the  interval  (0,1);  in  fact,  it  is 
n 

easily  shown  that 


/ 

o 

1 

/ 


1 

f  f  d£  ■  0 ,  /  f  "f  "dC  *  0,  m  i  n 
n  m  I  n  m 

0 


(3D 


The  fact  that  the  functions  f  (5)  satisfy  the 

n 

makes  then  ideal  approximating  functions  for  use  in 
Eqs.  (17)  and  (22).  Presumably,  a  solution  of  Eqs. 
approximated  by 


end  conditions  (Eq.  1 8 ) , 
numerical  attacks  on 
(17)  or  (22)  can  be 
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(32) 


yU.t) 


N 

E 

n*l 


qn(OfnU) 


Th.s  approach  reduces  the  problem  to  one  of  a  finite  number  of  degrees  of 
freedom.  The  generalized  coordinates  are  ^(t).  The  Lagrange  theory  of 
linear  vibrations  is  then  applicable.  Accordingly,  \et 


lO 


n 


(33) 


Then,  Eqs.  (24)  becomes 


♦  U) 


n  n 


(34) 


Since  L  depends  on  t.  note  that  ^  depends  on  t.  Consequently.  ^  is  not 

exactly  the  frequency  of  a  harmonic  motion.  The  general  solution  of  Eq. 
(34)  is 


gn(t*  V  0n> 


(35) 


in  which  An,  Bn  are  constants  of  integration.  An  infinite  series  solution 
of  Eq.  (22)  is  then 


y(5,t)  An,  (36) 

n-l 

Equation  (36)  automatical ly  satisfies  boundary  conditions  (Eqs.  IS)  at  the 

free  end  and  the  root.  One  might  expect  that  the  constants  (A  .  8  )  can  be 

n  n 

chosen  to  satisfy  the  conditions  at  t  *  0  and  t  -  t j ;  i.e.. 


10 


00 


E  V^(0'  V  6n>  •  F0(« 

n«  I 

09 

E  f„«W  W 

n«l 


(37) 


Al  ternat i  vel  y ,  the  constants  A^,  Bn  might  be  chosen  so  that  the  beam  has 
a  given  initial  deformation  (the  first  of  Egs.  37)  and  given  initial 
velocities;  i.e.,  y£(5,0)  might  be  specified  instead  of  yU.tj). 

If  L(t)  is  a  sufficiently  slowly  varying  function  of  t,  the  approxi¬ 
mate  solution  of  Eq.  (3*0  is 

g  “A  sinto  t  B  cosio  t  (38) 

*n  n  n  n  n 

Then,  Eq.  (36)  becomes 

00 

yU.t)  -  ^  (An  sinuint  ♦  Bp  cosu>nt)  f^U)  (39) 

n«l 


The  initial  value  problem  then  reduces  to 


Eb  f  U)  -  F  U) 

n  n  o 

n-l  (1*0) 

09 

E  *  ♦„«> 

n«l 

in  which  ^Q(^)  •  Yj(^*0)i  and  <^n°  i*  the  initial  value  of  Under  fairly 

broad  conditions,  the  constants  A  ,  B  can  be  chosen  to  satisfy  Eq.  |40) 

n  n 

when  F  (£)  and  i>  (£)  are  prescribed  functions.  Then  Eq .  (39'  represents  the 
o  o 

solution  of  the  initial  value  problem. 


I  1 


The  kinetic  energy  and  the  strain  energy  are,  respectively. 


2  2 
w  (A  cosu)  t  -  B  sinu)  t) 
n  n  n  n  n 


y5Ed<;  ■ 


El 

17 


8.  (A. 


s  i  nui  t 
n 


+  B 


COSU)  t)‘ 
n 


I 


Hence,  in  view  of  Eq.  (33),  the  total  mechanical  energy  is 


T  +  U 


n«l 


8  "(A/  +  B_  ) 


Thus,  the  total  mechanical  energy  is  independent  of  t,  provided  that  A 

n 

3/2 

and  vary  as  L  .  Apparently,  for  the  conditions  stated,  this  require¬ 
ment  ensures  conservation  of  mechanical  energy. 

(b)  F ( t )  Prescribed.  If  one  prescribes  F(t),  then  Eqs .  (17)  and  (20) 
subject  to  Eqs.  (18)  and  (21)  must  be  solved  simultaneously  in  order  to 
determine  L(t).  The  method  of  solution  for  both  cases  L ( t )  prescribed 
and  F(t)  prescribed  is  outlined  in  the  following  section. 


Method  of  Solution 

Equations  (17)  and  (20),  together  with  the  initial  and  boundary  condi¬ 
tions  given  by  Eqs.  (18)  and  (21)  form  a  system  of  two  nonlinear  partial 
differential  equations  in  the  two  unknowns  L(t)  and  y(x,t).  An  approximate 
numerical  solution  of  these  equations  can  be  obtained  by  a  Galerkin  finite 


element  formulation. 


Thus,  in  accordance  with  the  Galerkin  FEH,  we  form  the  approximate 


solutions  of  L  and  y  as  $  and  v  which  are, 


L( t)  -  4>(t)  ■  S ( t) 
y(x,t)  •  v(x,t)  -  NT(x)a(t) 


(41) 


where  N  are  a  set  of  quadratic  shape  functions  with  local  support.  A 
quadratic  shape  function  is  associated  with  each  nodal  point  of  the  dis¬ 
cretized  finite  element  model.  Quadratic  interpolation  over  an  element 
requires  each  element  has  three  nodal  points,  say  one  at  each  end  of  the 
element,  and  a  nodal  point  in  the  center  of  each  element. 

Figure  4  shows  that  the  odd  numbered  shape  functions  ,  ..., 

associated  with  nodal  points  at  the  ends  of  an  element  span  five  nodal 
points  (two  elements),  while  the  even  numbered  shape  functions  N^,  N^,  . 
span  the  three  nodal  points  of  one  element. 

In  accordance  with  the  Galerkin  FEM,  we  form  the  residual  function 
for  Eqs.  (17)  and  (20)  as. 


•  2  -2 

Rj (x,t)  »  v""  +  p[v"E2  +  2v'£ 


-  2v'£  1  -  v'C  f  +  v] 


(42) 


R2C,*tl  •  ♦  +  2?r  v"2(,’t>  -  Pt-  I  v'd.t)  +  v(!  ,t)}2] 

o  $ 


(43) 


The  finite  element  equations  are  obtained  by  requiring  that  the  residual 
functions,  given  by  Eqs.  (42)  and  (43),  be  orthogonal  to  each  of  the  basis 
functions.  That  is, 
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where  H  is  the  function  given  by 


e  e  e 


p*’  a2 


where  BT  denotes  a  boundary  term. 
Letting , 


,»•  8  4  . 2 

a  *  (7T  a2n-l  ‘  77  a2n  *  ~~2  ^n+l* 

le  i  ie 

c 


6  •  {f  (r  a2n-l  ■  r  a2n  *  a»-*i)  * 

e  e 


Equations  (*»9)  and  151)  become 


BT  ♦  Aa  -  Aa  ♦  2  Ca  -  2  T  c£  ♦  D? 

8  ~  8  ~ 


*  *  isr  7  l* ' ob]  *  sr FU) 

op  o 

Defining  the  following  terms, 
y  »  -(BT  ♦  Aa) 


•  —j  8  Aa  ♦  2  —j  Co 

8  8 

•2f c 


■Fc? 

ST  Fttl  -  ST  li  '  cb) 

o  o 


-  0 


(52) 


(53) 


(5M 


G  - 


where  8  is  the  value  of  8  at  the  previous  time,  Eqs.  (52)  and  (53)  become. 
Da  +■  e8  +  C*a  ♦  ■  y 

'  ~n  ...  (55) 

8  ■  S 

or,  in  matrix  form 

1  2n+l  2n+2 
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Equations  (55)  (or  (56)),  subject  to  initial  conditions  may  be 
solved  by  numerical  integration. 


Region  ft 

Differential  e*.  (6) 


ftlCM.ENft  Mint 


uhm.  MIT> 


ftNP-  CUM.  IM#»Ak  POINT 


GNP  “  GLOBAL  NOBAk  POINT 
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